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ABSTRACT 

We analyzed radio and X-ray observations of GRS 1915-1-105, between May 1995 and June 2006, 
focusing on the times characterized by radio flares and cycles of hard dips - soft spikes in the X-ray 
lightcurve. Assuming these flares to be discrete ejections, we applied a plasmon model to the radio 
data, with good agreement with the lightcurves. We fitted a total of 687 radio flares with a standard 
model of a plasmon. We found that the distribution of width is io = 1160 s with an rms deviation of 
360 s, while that of the amplitude is Smax = 59 mJy with an rms deviation of 28 mJy. The distribution 
of width is thus rather peaked, while that of the amplitude not. 

Regarding radio and X-ray links, this study confirms previous observations on smaller datasets, 
namely that X-ray cycles of hard dips - soft spikes are always followed by radio flares. A strong 
correlation is found between the length of X-ray "dips" in the X-ray lightcurves, and the amplitude 
and fluence of the subsequent radio oscillations. A model of an exponential rise of the form L15 GHz(At) 
= Lmax{^ — exp(— (Ai — /Strain) / t) is in good agreement with the observations, with the maximum 
fluence Lmax of the order 70 Jy.s, and the characteristic time r of the order 200-500 s. We discuss 
possible physical interpretations of this correlation, regarding the nature of the ejected material, and 
the physical process responsible for the ejection. 
Subject headings: X-rays: individual: GRS 1915+105 - X-rays: binaries - radio: jets 



1. INTRODUCTION 

K>" ' One of the main issues of Black-Hole Binary (BHBs) 
''nI" , physics is the link between accretion and ejection pro- 
^^ ■ cesses. Indeed, if the strong gravitational potential fa- 
f^ ' vors naturally the accretion of matter, the occurrence of 
^^ , relativistic ejections of plasma in these systems is still 
l/~j ■ far from being understood. A very promising object to 
O ' conduct this study is GRS 1915+105. 
f^ Indeed, GRS 1915+105 is a unique source in many re- 

spects. On the one hand, its main features classify it 
as a microquasar. GRS 1915+105 is a low-mass X-ray 
binary (LMXB), composed of a K- type star orbiting a 
V \ Black Hole (BH) of 14.0 ± 4.4Mq ()Harlaftis fc Greineii 
;_( ' |2004[) . As in most LMXBs, its X-ray spectrum can 
C^ ' be described with a soft disc blackbody (with kT ~ 
1-2 keV) , and a hard power l aw extending to >200 
keV (e.g. Rodr iguez et al.ll2008bl ). These are interpreted 
as an accretion disc and a Comptonizing region (of- 
ten called "corona"), respectively. At times, it displays 
superluminal ejections with true bulk velocities >0.9c 
(Mirabel & Rodriguez 199 4), or a ste ady compact jet 
(JDhawan et al.ll2000t iFuchs et al.ll2003[) . This overall de- 



scription corresponds to the "canonical" microquasar. 

On the other hand, GRS 1915+105 shows a 
unique wealth of beh avior (see, e.g., the review by 
iFender fc Bellonll I2004D . Contrary to transient micro- 
quasars, which are usually active for a few months be- 
tween quiescent states, GRS 1915+105 has been con- 
stantly active for the 17 years since its discovery with 
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GRAN AT (jCastro-Tirado et al.|[T99l . It shows a very 
high level of variability, which may be linked to the fact 
that GRS 1915+105 has a very hig h accretion ra t e com - 
pared with most other LMXBs. iBelloni et all (|20000 
classified its variability into 12 separate classes, which 
have been observed to recur almost identically over the 
years. In particular, the X-ray lightcurves display cy- 
cles of hard dips - soft spikes specific to GRS 1915+105. 
They are characterized by phases of low X-ray luminosity 
and hard spectra, which last between ^10 and ~2000s 
(see Fig. [1] and [5] for sample lightcurves) . These phases 
are ended by a short spike ('-^lOs), which marks the end 
of the hard phase. Note that because of the scaling, these 
spikes, although present, are not distinguishable on Fig. 
m right. Such features (hard X-ray dip ended by a spike) 
will be labeled simply "cycle" in the following. After the 
spike, the spectrum of the source becomes softer, and 
the X-ray luminosity decreases for a while. Usually, the 
light curve then starts to increase again and shows a sec- 
ond spike, longer (~100s), softer and brighter than the 
one marking the end of the hard dip. This second spike 
marks the beginning of a soft and highly variable phase, 
which ends when a new hard dip begins. 

The 12 GRS 1915+105 classes can, in turn, be inter- 
preted as transitions between three basic states, labeled 
A,B a nd C. Using the classification of iHoman fc Bellonil 
(|2005D . these 3 states can be related to the canonical 
Soft State, Soft Intermediate State and Hard Intermedi- 
ate State, respectively. Cycles can thus be described as 
short C states separated by sequences of A and B states. 
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Fig. 1.— Simultaneous X-ray and Rylc observations of GRS 1915+105. On MJD 51863.64, GRS 1915+105 was in the variability class v 
(left), on MJD 53324.7 it was in class A (right). The plots show on the top the 2-60 keV RXTE lightcurve (panel (a)), and the 3-30 keV 
JEM-X lightcurve (panel (d)). Panels (b) and (e) display the 15 GHz radio lightcurve, while panels (c) and (f) display the hardness ratios. 
The vertical dotted lines mark the beginning and the end of the X-ray dips on each panel. The horizontal dotted lines are the estimated 
level of background radio emission. On top of these lines, the continuous lines show the fit of the plasmon model on the radio data. Both 
observations correspond to a confidence index of 1 (see section r4.3l l. 
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Fig. 2. — Simultaneous RXTE and Ryle observations of GRS 1915+105, during variability class /3. Panels and lines are identical to Fig. 
[l] On the l eft side, the radio lightcurve corresponds to a confidence index of 2, while on the right side it corresponds to an index of 3 (see 
section l4.3| l. 
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In the radio range, GRS 1915+105 displays also a very 
high level of variability, which is related to the 12 vari- 
ability classes. In particular, strong radio oscillations 
or isolated flares a re detected, which are clearly related 
to the X-ray dips (jPoolev fc Fendedll997l: IMirabel et all 
[1998) . The power required to generate the synchrotron 
emission of these oscillations is likely to be a very signif- 
icant fraction of the entire accretion ene rgy of the sys- 
tem (jFender fc PoolevI [2000l: iMeieii [200l . Radio oscil- 
lations arc particularly strong during /3 and v classes 
(JKlcin-Woh ct al. 2002). The a, A, 6*, n and p classes 
are characterized by a weak or variable radio flux. 

Since X-ray cycles and radio oscillations seem to occur 
at the same times, it is temptin g to look for a link be- 
tween them. Previous analyses (jPoolev fc Fended 119971 : 
iKlein-Wolt et al.l 120021 ) have indeed found such a link, 
namely that X-ray cycl es are followed by radio flares dur- 
ing /3 and v classes. iRodriguez et al.l f|2008a) showed that 
this behavior can be extended to A classes, and that the 
time lag between the X-ray spike at the end of a cycle 
and the maximum of the subsequent radio flare is roughly 
constant, regardless of the characteristics of the preced- 
ing dip. This led to the interpretation that the spike is 
the trigger of the radio flare. 

This possible X-ray - radio link, which is specific to 
GRS 1915+105, can be related to the behavior of "con- 
ventional" transients. Indeed, the spike also marks an 
abrupt transition from harder to softer spectral states, 
and is followed by radio activity. In the temporal domain, 
timing analyses suggest that this transition i s associated 
to ty pe B quasi-periodic oscillations (QPOs. lSoleri et all 
|2008[) . These sequence of events and timing behavior are 
in analogy with the detection of major radio flares and 
type B QPOs at the time of X-ray spect ral transitions in 
many microquasars ([Fender et al.ll2009t ). 

The purpose of this paper is to explore the possible link 
between X-ray dips and radio flares, with an extended 
data set and a physical treatment of the radio ejections. 
Our analysis is centered on a, /3, A and v classes, which 
are well covered by observations. A few words will be 
said about the other cyclic classes in section l475l 

Section 2 starts with a description of the available 
data, and the reduction process. Section 3 describes the 
method and models used to analyze the radio data, while 
section 4 explores the links between radio and X-ray be- 
havior. In section 5, we discuss the physical interpreta- 
tion of these links. 

2. OBSERVATIONS AND DATA REDUCTION 

In the radio range, the Ryle Telescope (RT) provides an 
extensive coverage of GRS 1915+105 between May 1995 
and June 2006. Its sensitivity (a few mJy with a 32s 
temporal binning) is sufficient to follow discrete ejection 
events in GRS 1915+105 and adjust physical models to 
represent them. 

2.1. Radio data: Ryle Telescope 

RT obse rvations follow the scheme described by 
iPoolev fcl^ender (1997) . The setup measured the Stokes 
parameter I + Q, calibrated using the phase calibrator 
B 1920+154, and the amplitude calibrators 3C 48 or 
3C 286. The resulting 15.2 GHz light-curves were binned 
at 32s. Using data taken during a quiet period, we es- 
timated a standard deviation of 5.7 mJy, and used this 



value as nominal error bars. 

2.2. X-ray data: RXTE/PCA and INTEGRAL/JEM-X 

The RXTE data were reduced wi th the HEASOFT 
v6.5 so ftware package, as described in [Rodriguez et al.l 
([2008al) . We extracted Is resolution light curves from 
the Proportionnal Counter Array (PCA) in the 2-60 keV 
range, as we ll as in the two energy bands defined in 
IBelloni et all HOOO): 2-5.7 keV and 5.7-14.8 keV. We 
then computed the Hardness Ratio HR=5. 7-14. 8/2-5. 7 
keV. 

We also used data from t he INTEGRAL monitor - 
ing campaign described in IRodriguez et al.l ([2008a[) . 
The INTEGRAL data were reduced using the standard 
Off-line Scientific Analysis (DSA) v7.0 software 
package ([Goldwurm ct al. 2003). Data from the JEM- 
X instruments were reduced to construct 10s lightcurves 
in the 3-7, 7-15 and 3-30 keV energy range (channels 
46-95, 95-159 and 46-210 respectively). Since JEM-X is 
less sensitive than the PCA, the Hardness Ratio (HR) is 
too noisy to be more than a hint to the classification of 
variability classes. 

3. A LOOK AT THE RADIO FLARES 

In the radio range, GRS 1915+105 displays a very 
rich behavior, characterized by the presence of discrete 
events. Two types of flares are visible. Firstly, a few gi- 
ant flares have been detected, which last several days and 
reach hundreds of mJy (Fig[3|). Then, numerous flares, 
with a typical amplitude of a few tens of mJy are also vis- 
ible, which last around 1 hour. These smaller flares are 
either isolated or repetitive, and several successive radio 
flares can overlap in the radio lightcurve (Fig. [2|). Before 
looking at X-ray observations and variability classes, we 
tried to characterized these small flares, i.e. their width, 
amplitude and recurrence period. 

To do this, we searched for radio oscillations in the 11 
years of RT data. They contain about ~1000 clearly- 
identifled flares. In order to study these flares, we need 
an emission model to map their evolution. 

3.1. Flare models 

Multiwavelength studies provided good confidence that 
synchrotron processes are responsible for the emission 
during radio flares. However, doubts remain on the exact 
nature of the synchrotron emitting medium, and on its 
geometry. Th e most widely used model was described by 
Ivan der LaanI ([l966i) . who calculated synchrotron emis- 
sion from an adiabatically expanding cloud. This model 
was first developed for Active Galactic Nuclei (AGNs), 
for which it is still commonly used. In the case of 
GRS 1915+105, IMirabel et all p998[ ) observed delays 
in the lightcurves at different radio and infrared wave- 
lengths compatible with this plasmon model. 

Another possibi lity was described by 

iHiellming fc JohnstonI (|l988[ ). Starting with the 

van der Laan model, they calculated the emission from 
conical twin jets, and adequately predicted the emission 
of several X-ray binaries. They also provided a geomet- 
rical refinement for the van der Laan model, leading 
to slightly different radio lightcurves for expanding 
s pherical bubb l es. 

[Kaiser et al.l ([2OOO0 described an internal shock model 
for the observed synchrotron emission. In this scheme, a 
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Fig. 3.— Radio lightcurve of the MJD 50749-50759 major flare 
seen by the Ryle Telescope. The insert shows a zoom on the 
hghtcurve, where smaller flares are clearly visible on top of the 
major event. 

quasi-continuous jet emission is "lit-up" by shock fronts 
traveling along the jets. This model led to similar results 
as the plasmon model, although with a smaller power 
emitted by the central engine. 

Up to now, no scenario is clearly favored by the ob- 
servations. In this paper, we are first interested in the 
morphological parameters of the flares, i.e. amplitude, 
width and fluence. For these parameters, every model 
will produce approximately the same values, as long as 
it describes accurately the observations. Therefore, we 
used the simplest one in our study, i.e. the van der Laan 
plasmon model, but it should be noted that the physi- 
cal nature of the flares could be different from a simple 
spherical plasmon. 

3.2. Description of the model adopted 

In the Ivan der LaanI (|1966f l model, flaring emission is 
produced through the adiabatic expansion of an initially 
optically thick blob of synchrotron-emitting relativistic 
electrons. During the first phase, the increase in the 
blob's surface area causes an increase in the flux. Then, 
the curve turns over as the plasma becomes optically thin 
because of the reduction in magnetic field, the adiabatic 
cooling of the electrons, and the reduced column den- 
sity as the blob expands. Simultaneous observation at a 
lower frequency will show the same fractional rate of in- 
crease, but the maximum will be reached later and have 
a smaller value. 

As a starting point, we assume that the electrons have 
an energy spectrum n[E) oc E~^ . The synchrotron opti- 
cal depth at frequency v of an homogeneous and spherical 
blob then scales as 



Ty{t) = To 



-(p+4)/2 



m 



-(2p+3) 



where R{t) is the radius of the spherical blob at a given 
time. The flux density then scales as 



Siy{R) — So 



5/2 



R 

Rq 



1 — exp(— t) 
1 - exp(-ro) 



^u{i) — ^max 1 , I 1 t 



Here, i?o, 5*0 and tq are the size, flux density and 
optical depth at the peak frequency of the synchrotron 



spectrum vq (|van der Laan|[T966f) . tq only depends on p 
through the condition 

e^° -To(p + 4)/5-l = 

Thus, given the particle energy spectral index p and 
the peak flux 5*0 at a given frequency vq, this model 
predicts the variation in flux density at any other fre- 
quency as a function of the expansion factor {R/ Rq). A 
model for R{t) is needed to map this relation to time: we 
will assume a simple linear expansion at constant speed 
Vexp- The major ejections detected from G RS 1915-1-105 
indicate an expanding velocity of ^ 0.8c (jFender et al.l 
|1999[ ) This model predicts a sharp increase of luminos- 
ity, followed by an exponentially decreasing tail. The last 
parameter, the energy spectral index p depends on the 
source at hand, and is usually taken between 1 and 3. 

In order to provide precise measurements of the plas- 
mon parameters, the model requires simultaneous mon- 
itoring at several frequencies. In particular, the peak 
frequency vq cannot be determined without a spectral 
monitoring. Since we observe at a fixed frequency of 
v = 15.2 GHz only, we will use a simplified version of the 
formula. Wc define the characteristic length of a flare to 
be 

D 1 /,, ^ (P + ^). 
, -n-O (2p + 3) / '^O \ 145+61 
to - Tg I 1 

Vexp ^V ^ 

This leads to the final fitting function 

t " 

to 

Thus, the model will have 3 free parameters: the char- 
acteristic length of the flare tg , the maximum amplitude 
Smax and the initial time of the ejection. 

Several radio observations show a more complex be- 
havior: distinct radio oscillations are visible on top of a 
strong, slowly varying "background flux" . Fig. |3] shows 
an example of such observations: a very bright flare oc- 
curs between MJD 50749 and 50759, on top of which 
smaller flares are clearly distinguishable. In order to iso- 
late the contribution from individual flares, one thus has 
to remove the contribution from the background flux. 
To do this, we adopted an empirical solution: since the 
additional flux varies on a timescale larger than a few 
radio oscillations (see Fig. [5]), we assumed the back- 
ground flux to be linear over a few flares. Thus, we re- 
moved from the radio flux a linear function of the form 
Sback{t) = At + Kq, where A is computed by a linear 
fit to the radio data, and Kq is let free under the con- 
dition that Sback remains below the minima of the radio 
lightcurve. Under this definition, Sbackit) represents the 
source flux underlying the small flares: a fraction of the 
flare flux will be incorporated in this quantity, which will 
lead to a slight underestimation of the flare fluxes. 

3.3. Typical radio flare parameters 

Once applied to the observations, the plasmon model 
adequately fits all data, with a reduced chi-squared xt 
between 1 and 3. This is especially true for isolated flares, 
when no strong background correction is needed. Fig. [U 
panels (b) and (e), shows such fits, with xt = 1-26 and 
1.79 for the two flares peaking on MJD 51863.68 and 
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Fig. 4. — Distribution of the characteristic width and amplitude 
of the radio flares. The histograms are computed using 687 flares. 
The continuous lines are gaussian fits to the histograms. 

MJD 53324.79, respectively. These good results give us 
good confidence in our approach. 

We checked for the influence of the spectral energy in- 
dex p of the electrons for the predicted lightcurve. p = 2 
and p = 3 gives very similar results: the shape of the 
predicted flare is similar, and adequately fits the obser- 
vations. Using p = 1, however, the predicted flare is 
narrower, and fails to reproduce the observations. In 
the following, the results are given for the nominal value 
p = 2. 

In order to extract the main characteristics of the 
flares, we applied the plasmon model systematically to 
the RT data. In doing so, we excluded very noisy data 
or incomplete flares. A total of 687 flares have been fit- 
ted. Fig. m presents the distribution of their characteris- 
tic width (given by the to parameter) and amplitude, to 
varies between and 3000 s, while the maximum ampli- 
tude varies between 10 and 150 mJy 

The distribution of flare width and amplitude appear 
to be gaussian. Gaussian fits to these histograms give 
to — 1160 s with an rms deviation of 360 s, and a mean 
amplitude of Smax — 59 mJy with an rms deviation of 
28 mJy 

Note that the distribution of these parameters is bi- 
ased by four main uncertainties. Firstly, very small 
flares will hardly be distinguishable from the background 
noise, which varies between 5 and 6 mJy on the 32s bin 
lightcurve. Thus, the histograms are truncated at small 
values. Then, the uncertainty in the background correc- 
tion will tend to increase the width of the distribution: 



the linear model is a strong approximation. 

In the case of successive flares, an overlap between 
them can occur. This was first not taken into account in 
the background correction. The uncertainty this intro- 
duces will tend to shift the maximum of the to and Smax 
distributions towards higher values. The effects of over- 
lapping flares are considered in the next section. Finally, 
in the few cases where, because of the simplified nature 
of the model, the fit leads to a rather high xt around 
3, the fiare parameters are not well constrained. In this 
case, the effect on the distributions is an increase of their 
widths (rms deviations). 

4. RADIO - X-RAY LINKS 

4.1. Relevant observations 

We cross-correlated the log of the radio pointings to X- 
ray observatories, and found a total of 352 RXTE point- 
ings and 7 INTEGRAL observations simultaneous with 
RT data. 

Since our purpose is to study the radio behavior during 
X-ray cycles, we again filtered our dataset to retain only 
those observations with good simultaneous coverage. We 
looked for observations with (1) enough X-ray coverage 
to detect clearly a hard X-ray dip and (2) radio coverage 
during the following hour to detect a possible activity. A 
total of 28 RXTE pointings and 3 JEM-X observations 
display these characteristics for a, 13, A and i' classes. 
These pointings contain 54 hard X-ray dips with radio 
coverage, which are listed in table [TJ 

Fig. [U and H display four examples of RXTE /PCA 
and INTEGRAL/ JEM-X hghtcurves with X-ray dips, for 
classes /3, v and A. The corresponding HRs displayed in 
panels (c) and (f) illustrate the sharp spectral hardening 
that occurs during dips. The end of each dip is marked 
by a short spike, particularly visible in /? and v classes. 
In the radio range, oscillations or isolated fiares follow 
this spike. 

4.2. Consecutive flares 

In order to compare the flare characteristics to the X- 
ray data, our simple model needed to be refined, no- 
tably in the case of consecutive flares. Indeed, several 
radio lightcurves of GRS 1915-1-105 are characterized by 
quasi-sinusoidal modulations, which can be interpreted 
as repeated discrete ejections. In this case, the emission 
coming from different flares overlaps in the lightcurve. 
In order to disentangle the emission from each one of 
them, we modeled these successive ejections using multi- 
ple identical flares, separated by a constant interval tiag- 
We checked that, when looking at a sequence of a few 
radio flares, the time lag between two maxima remains 
constant within a few % (Fig. ^. 

This "multiflare" approach can be compared to the 
previous approach used on the 687 flares. By taking the 
overlaps into account, the refined model leads to values 
that are lower by at most 10 % for the amplitude, width 
and fluence of a given fiares. This has to be compared to 
the uncertainties on these parameters, which lie between 
~10-50 %. Therefore, in the general study above, this 
effect is not crucial in the determination of the global 
characteristics of the flares. Since the refined model also 
leads to lower xt values, the main source of uncertainty 
that now remains is the level of the background radio 
fiux. 
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TABLE 1 
List of observations with simultaneous X-ray and radio 

COVERAGE, during CLASSES a, H, X AND I/. 



Start Date 


Class 


X-ray 


Dip length 


Confidence 


(MJD) 




instrument 


(s) 


Index 


5038L495 


u 


PCA 


1526±3 


2 




V 


PCA 


1333±5 


2 




V 


PCA 


> 1800 


2 


50698.658 


/3 


PCA 


562±5 


1 




/3 


PCA 


554±1 


1 




/3 


PCA 


542±1 


1 


50751.687 


/3 


PCA 


767±1 


3 




/3 


PCA 


> 515 


3 


50751.754 


/3 


PCA 


766±1 


3 


51343.042 


/3 


PCA 


> 410 


3 




/3 


PCA 


1000±2 


3 




/3 


PCA 


> 440 


3 




/3 


PCA 


884±2 


3 


51352.957 


/3 


PCA 


> 1430 


- 




/3 


PCA 


1411±2 


- 


51353.025 


/3 


PCA 


> 1400 


3 


51476.639 


/3 


PCA 


756±3 


3 




/3 


PCA 


> 400 


3 


51814.727 


V 


PCA 


> 1830 


3 


51863.594 


V 


PCA 


1125±1 


1 


51863.661 


V 


PCA 


1009±1 


1 


51870.554 


p 


PCA 


> 540 


2 




p 


PCA 


1076±2 


2 


51870.627 


n 


PCA 


948±2 


3 




p 


PCA 


> 530 


3 




n 


PCA 


955±2 


3 




p 


PCA 


> 740 


3 


51877.522 


V 


PCA 


1470±10 


- 




V 


PCA 


> 700 


- 


51877.599 


V 


PCA 


1828±10 


2 




V 


PCA 


> 1400 


2 


51885.509 


n 


PCA 


925±75 


3 


51954.302 


V 


PCA 


1140±5 


2 


51954.441 


V 


PCA 


> 815 


1 


52108.806 


V 


PCA 


> 1800 


2 


52108.877 


V 


PCA 


> 1400 


2 


52171.799 


A 


PCA 


140±2 


1 




A 


PCA 


72±2 


1 




A 


PCA 


82±2 


1 


52500.893 


a 


PCA 


1975±20 


- 


53150.057 


P 


PCA 


857±3 


3 


53296.373 


V 


JEM-X 


1826±20 


1 




V 


JEM-X 


1727±20 


1 


53296.728 


V 


PCA 


> 1400 


1 


53296.794 


V 


PCA 


> 1000 


1 


53318.559 


p 


PCA 


> 200 


1 




p 


PCA 


794±2 


1 


53324.277 


A 


JEM-X 


130±20 


1 




A 


JEM-X 


180±10 


1 




A 


JEM-X 


600±20 


1 


53503.633 


P 


JEM-X 


420±20 


2 




P 


JEM-X 


460±10 


2 


53703.573 


V 


PCA 


> 700 


2 


53703.639 


V 


PCA 


827±10 


2 



4.3. Confidence indices 

However, in the case of a strong background, this level 
is poorly constrained. Indeed, there is a degeneracy be- 
tween the offset level Kq and the other parameters: if Kq 
is fixed to lower values, the fitting routine will change the 
other parameters to increase the overlap between fiares, 
which will lead to an acceptable fit. This effect is partic- 
ularly true in the case of strong overlap. 

In order to quantify this issue, we attributed a "confi- 
dence index" to our data. This index is based on two cri- 
teria: the level of background flux, and the morphology 
of the observed flares. Indeed, when a radio oscillation is 




400 600 SCO 1000 1200 1400 
Lag (s) 

Fig. 5. — Distribution of the delay between the X-ray spike and 
the maximum of the subsequent radio flares. The histogram has 
been computed for the 34 flares for which the spike was visible. 
The continuous line is a gaussian fit to the histogram. 

fully shaped, in particular with a clear exponential tail, 
we have good confidence that the background subtraction 
is accurate. On the other hand, in the case when only 
the "tip" of the flare is visible, the background level can- 
not be precisely constrained, which is a source of higher 
uncertainty for the determination of flare properties. 

We attributed a confidence index of 1 to the most re- 
liable data: observations with fully shaped fiares, either 
isolated or on top of a low background emission (<15 
mJy, Fig. [T]). An index of 2 corresponds to observa- 
tions where the exponential decrease is still clearly visi- 
ble, with moderate background emission (<30 mJy, Fig. 
[21 left). Finally, observations with almost sinusoidal os- 
cillations on top of a strong background (up to 120 mJy, 
Fig. [21 right) correspond to index 3. Thus, parameters 
deduced from index 3 observations should be considered 
with care, while index 1 observations would produce re- 
liable parameters. 

Among the observations listed in Table [H four obser- 
vations had to be excluded from our data set: on three 
observations, two successive flares are visible but not dis- 
tinguishable, while on one more observation the radio 
data are too noisy to get any constraint on the fiare. 

4.4. X-ray dips and occurrence of radio flares 

To measure the characteristics of the X-ray dips, we 
used the following definition. Time is defined as the 
time when the phase of highly variable X-ray fiux ends. 
This time is also the first point of spectral hardening 
(visible in the HR). We used the position of the maximum 
flux of the X-ray spike to mark the end of the dip, whose 
duration will be noted At. At this time, the HR has gone 
down to values close to the prc-dip phase. 

In order to measure the amplitude of the final spike, the 
difficulty lies in the determination of the beginning of the 
spike: there is no particular feature marking this point. 
Instead, we used the minimum flux of the dip, whose 
value is close to the fiux near the end of the dip. Thus, 
the amplitude was taken to be the difference between the 
top of the spike and the minimum flux of the dip. 

Using the parameters given by the fits, we looked 
for links between the X-ray and radio behaviors of 
GRS 1915-1-105. We confirm the association between X- 
ray dips and radio flares: among the 54 X-ray dips listed 
on table 1, 52 are directly followed by a radio flare, only 
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Fig. 6. — Width of the radio flares (to), as a function of the 
length of the preceding Hard X-ray dip At. Blue and black circles 
correspond to confidence index f observations, red squares to index 
2, and green triangles to index 3. The error bars are estimated 
at the 3(T confidence level (note that errors on index 3 data are 
questionable) . 

the 2 shortest dips being not fohowed by a detectable 
activity. We thus confirm that X-ray dips during a, /3, 
A and v classes are always followed by radio flares. Fig. 
[5] shows the distribution of the delay between the X-ray 
spike and the peak of the radio flare, for 34 X-ray dips. 
For 18 other dips, a flare is also detected but the date 
of the spike is not known because of the lack of X-ray 
coverage. 

The average delay between the X-ray spike and the 
peak of the radio flare is 1040 ± 185 s (at Ig), a value 
consist ent with the one calculated by [Rodriguez et al.l 
(j2008aD . As for the derived start time of the ejection, 
according to the van der Laan model, the ejection of 
matter is coincident with the time of the X-ray spike 
within less than 5 minutes. 

4.5. Other cyclic observations 

Since this connexion between X-ray dips and radio 
flares seems ubiquitous, it is interesting to look for cy- 
cles not followed by detectable flares. Firstly, very 
short cycles, such as the two A cycles observed on MJD 
52171.799, are not followed by a detectable radio activ- 
ity (upper limit of ^2 mJy) . Given that A cycles longer 
that 100 s, such as the one displayed on Fig. [U right, 
are followed by a weak radio flare, an explanation may 
simply be that such flares trigger too small ejections to 
be detected by the Ryle Telescope. 

Then, the 0, k and p classes are not followed by distinct 
radio flares, either. These classes were not included in the 
analysis so far. k and p classes display very short cycles, 
between 10 and 50 s. These cycles recur very quickly, 
with less than 100s between two consecutive X-ray spikes. 
Our fits show that the radio fiares of GRS 1915-f 105 peak 
around 1000 s after the end of the cycle, and last typically 
1200 s, thus they could not be seen individually during 
these classes. Note that a weak radio fiux is detec table, at 
~3-4 mJy (as already seen bv iKlein- Wolt et al.l ([2002)), 
which could be compatible with the sum of very small 
radio ejections occurring after each cycle. 

6 classes display longer cycles, with X-ray dips lasting 
typically between 300 and 600 s. These cycles also recur 
quickly, with less than 1000s between two consecutive X- 
ray spikes. In this case too, individual flares will not be 
distinguishable, but should produce a varying radio flux. 



Fig. 7. — Amplitude of the flare (normalization Smax), as a func- 
tion of the dip-length At. The symbols are identical to Fig. [S] 

Radio observations duri ng this class show a str ong ra- 
dio flux, at -10-80 mJy ()Klein-Wolt et al.l 120021 ). This 
flux is highly variable, but displays no specific pattern. 
This value of ~ 10-80 mJy is significantly lower than the 
maximum fluxes observed after /3 and v cycles of com- 
parable duration. Thus, this radio activity could be ex- 
plained by radio flares, but of lower amplitude. Note that 
9 cycles are also characterized by less pronounced X-ray 
dips: the Hardness Ratio is lower and the minimum flux 
higher than during other cyclic classes. This may explain 
their different behavior regarding radio activity. 

To put it in a nutshell, observations of 6, k and p classes 
are not entirely conclusive. On the one hand, they can be 
understood within the framework of radio ejections; the 
distinct emissions of matter are simply too quick to be 
separated in the radio lightcurve. On the other hand, the 
observed radio flux is lower than expected if we followed 
the trend given by /3, A and v classes. This probably 
means that several features determine the characteristics 
of the radio activity: the length of the X-ray dips may 
not be the only parameter required to explain all the 
radio flares. 

4.6. A correlation between duration of X-ray dips and 
fluence/ amplitude ? 

Now, let us look for correlations between the various 
parameters at hand. The amplitude of the final X-ray 
spike is not related to the width or amplitude of the fol- 
lowing radio flare: the maximum X-ray flux seems to be 
random. No link was found between the minimum flux 
during a dip and the characteristics of the following flare 
either. A slight link, although not statistically signif- 
icant, was found between the minimum X-ray flux and 
the duration of the dip: lower minimum fluxes tend to be 
rel ated to longe r dips. This link was already mentioned 
bv lBelloni et al.l (|1997| ) in the case of class k observations. 

Finally, a more pronounced trend is visible between 
the duration of the dips Ai and the parameters of the 
subsequent flares: longer dips seem to be followed by 
more important radio ejections (Fig. |6l [7] and [8]) . In 
particular, a correlation is visible between Ai and the 
amplitude Smax of the following radio flare. A correlation 
is also visible between Af and the fluence Li^ghz of the 
following radio flare. 

To quantify this, we used the Pearson and Spearman 
tests of correlation. The Spearman test detects any 
monotonic correlation in the data, while the Pearson test 
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TABLE 2 

Significance of the correlation between the lencjth of 

X-RAY DIPS Ai AND THE AMPLITUDE Smax, WIDTH io AND 

FLUENCE Lic,CHz OF SUBSEQUENT RADIO FLARES. VALUES ARE 

EXPRESSED IN LEVELS OF SIGNIFICANCE FOR, THE TWO TESTS. 



Quantities 


Test 


Index 1 


Index 1 & 2 


All data 


LlSGHz vs. At 


Spearman 


>99.99% 


98.4% 


98.3% 




Pearson 


87.3% 


65.1% 


56.5% 


Smax VS. At 


Spearman 


99.98% 


>99.99% 


99.92% 




Pearson 


85.6% 


80.5% 


59.0% 


to VS. At 


Spearman 


77.0% 


1.0% 


1.0% 




Pearson 


39.1% 


1.2% 


2.3% 




1000 

Dip length (s) 



detects linear correlation. Values of these indices are re- 
ported in Table [2] in confidence levels. The two tests 
show a strong correlation between the length of the X- 
ray dip At and the maximum of the following flare Smax ■ 
The correlation between the length of the dips and the 
fluence Li^cHz is also high. Spearman and Pearson tests 
are not significant in the case of length of the X-ray dip 
At versus width tp of the flares. 

Note, however, that the Li^qhz vs. At and Smax vs. 
At correlations are less pronounced when ignoring the 
five points corresponding to At < 200s. In the case of 
Smax vs. At, the level of significance for the Spearman 
test drops to values between 60% and 98%, for the "all 
data" and "Index I" subsets, respectively. Indeed, as 
visible on Fig. [71 the index 1 data are still correlated 
when ignoring the A points. In the case of Lisc-ffz vs. At, 
however, the non-A points are almost compatible with a 
constant. 

During the A observations, there is no radio emission 
prior to the X-ray spike, so the characteristics of the flare 
are well constrained. We therefore have good confidence 
on the robustness of these values. However, a crucial test 
for the reality of the Smax vs. At and Li^qhz vs. At 
correlations would be to observe X-ray dips in the range 
200-500 s. More observations are needed to populate 
this range in order to refine the link between X-ray dips 
and radio flares. 

Looking at the fluence vs. dip-length relationship, the 
correlation does not seem to be linear nor affine. Indeed, 
short X-ray dips produce much weaker radio flares than 
long dip s (see Fig. El). Th is is co nsistent with the state- 
ment of iKlein-Wolt et al.1 (|2002[ ) , who mentioned that 
dips shorter than lOOs do not seem to be followed by 
radio flares. 

Therefore, we looked for other functions to describe the 
relationship between the (radio) fluence L15 ghz and the 
length At of the X-ray dip than a simple linear or affine 
law (Fig. [8|). Two types of functions produced good 
results, which are listed in Table [3l Firstly, the data are 
well described by a function of the form ii5 GHz(At) = 
yl(At— Ato)'', with 7 between 0.17 and 0.35 depending on 
the subset considered. Then, an exponential function of 

the form L15 GHz(At) = Lmax{^ - e ' ) provided 

very similar x^, with a characteristic time r between 100 
and 500 s, and a saturating fluence Lmax between 55 and 
80 Jy.s. 

Both models are truncated at short dip-length, as vis- 
ible on the plots. In these two cases, the quantities Atg 
and Atmin correspond to the length of the X-ray dip 
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Fig. 8. — Fluence of the radio flares as a function of the dip-length 
At. The symbols are identical to Fig. [6] The top panel displays all 
available observations, while the lower one displays index 1 data 
only. The continuous line is an exponential fit and the dashed line 
a power-law fit to the data on each panel. 

below which no flare is detected (L15 ghz = 0). They 
measure Ato and Atmin between 50 and 130 s. This is 
compatible with the empirical observation that no flare 
are detected after dips shorter than ~100 s. 

5. DISCUSSION 

Our analysis extends t he previous result s on th e oc- 
currence of radio flares. IKlein-Wolt et al.] ([20 02) sug- 
gested that, during (3 and v classes, long and spectrally 
hard State C int e rvals a re followed by radio oscillations. 
[Rodriguez et alj ()2008aD extended this connexion to A 
classes, and identified the X-ray spike at the end of the 
X-ray dips as the trigger to the ejection. With an ex- 
tended data set, we confirm that, during a, /3, A and v 
classes, every X-ray dip in the lightcurve is followed by 
a radio flare. Moreover, the beginning of this flare coin- 
cides with the X-ray spike within less than 300s. On k 
and p classes, we lack the sensitivity to confirm or reject 
a similar behavior. Finally, 9 classes may also follow this 
scheme with, however, radio flares significantly smaller 
than during /3 and v dips of comparable duration. 

Using the available /3, A and v observations, a trend is 
visible between the length of the X-ray dips and the char- 
acteristics of the following flare: the longer the X-ray dip 
is, the bigger the following radio fiare will be. This result 
does not directly depend on the plasmon model used to 
characterize the radio data. Indeed, this model we chose 
was used mainly to distinguish between the energy emit- 
ted during each individual fiare in a sequence of radio 
fiares, rather than as a precise description of each fiare. 
Therefore, every "realistic" function with a fast rise and 
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TABLE 3 

Values obtained by fitting power-law and exponential functions to the dip 
length vs fluence data. the 2 functions used were of the form 

Af-At„ 
-^15 GHz{^*) = Lmaxi^ 



') AND Lis GHz(Ai) = A{At - Ato)'' . Errors are 



GIVEN AT THE 68% CONFIDENCE LEVEL (Icr). 



Set 


^max 

(Jy-s) 


(s) 


(s) 


xl 


A 
(Jy-s) 


(s) 


7 


xl 


Index 1 

Index 1 & 2 

All data 


80±13 
55±20 
63±10 


490±190 
90±67 
220±75 


50±70 

120±20 

90±40 


0.34 
2.09 
3.90 


6±4 

77±6 

9±2 


125±15 

130±2 

129±6 


0.35±0.10 
0.17±0.06 
0.28±0.02 


0.29 
2.11 
4.08 



slow decay able to fit the shape of the radio flares would 
lead to the same result. In particular, emission from a 
conical jet or a shock-heated compact jet would have al- 
most the same temporal shape, and thus lead to similar 
results concerning the basic parameters explored here. 

This link seems to be related to the length of hard X- 
ray dips: longer dips lead to more energetic ejections. 
Although the data are scarce, this link does not seem 
to be linear, and can be described, for instance, by a 
power-law or exponential relationship. Given the high 
uncertainties on the data, it is not possible to discrimi- 
nate between these functions. 

Yet, in the second case, one could think of an at- 
tractive scenario to explain the relationship. Obser- 
vations show that the beginning of the dip is marked 
by a quick increase in the accretion disk inner ra- 
dius. Then during the dip itself, the disk draws closer 
to the BH dBelloni et a l.|[l997l: iMig liari fc Bellom| [200l 
[Rodriguez et al.l l2008bD . This evolution brings more 
gravitational energy close to the BH. Let us suppose that, 
during the Hard X-ray dip, energy is somehow extracted 
from the accretion disk, and accumulated in the sur- 
rounding medium. Thus, the luminosity of the disk and 
corona slowly increase. At the same time, the amount 
of energy lost by the corona in a given time grows as 
the energy density of the corona increases. Thus, the to- 
tal energy stored inside the corona saturates after, say, 
^ 1000s. Then, at some point, this energy is released in 
the form of a quickly expanding blob of matter. In the 
X-ray ligh tcurve, this ejection is marked by a short X - 
ray spike (jMirabel et all 119981: iRodriguez et al.ll2008aD . 
Note that the ejected material can come from the corona 
itself, or from the inner accreting disk; either way, after 
the spike the corona is not visible anymore, and the disk 
is closer to the BH. 

In the case of a jet (either steady or shock- heated), 
the scenario is similar: during the dips an injection 
of material into the je t from the corona takes place 
(jKlein-Wolt et all 120021 ). Assuming a more or less con- 
stant injection rate, longer X-ray dips will thus be radio- 
brighter, as more material is ejected. Furthermore, dur- 
ing the dip, as the disk moves back in towards the com- 
pact object, one can assume that it reduces the injection 
rate, which would explain the apparent saturation in the 
flare fluence. Then, at some point, the disk is so close 
to the BH that it suppresses the injection of particles. 
The X-ray dip ends, and radio emission is detected af- 
terwards. 

However, the jet interpretation suffers from an impor- 
tant caveat: it predicts lightcurves different from the ob- 



served radio flares. In deed, using the conical jet from 
iHjellming fc JohnstonI (|1988i ) . it is possible to model the 
radio emission from a transient jet. During a given dip, 
the jet would be active for ^30 min. With a velocity of 
about Vexp — 0.8c, the jet would have a total extension 
of ~ 4.10^'^cm, or about ~ 10^i?s- The resulting syn- 
chrotron emission would thus emanate from electrons at 
very different temperatures, due to adiabatic expansion. 
When looking at a fixed radio frequency, the lightcurve 
would then be broader than that of a single plasmon. 
Besides, the radio emission would be detected from the 
beginning of the X-ray dip, or with a constant time-lag 
with it. Since the data show that the radio emission be- 
gins within a few minutes from the X-ray spike at the end 
of the dip, with no dependance on the length of the dip, 
this interpretation seems less probable than the discrete 
ejection model. 

It is interesting to note that the distribution of width 
of the radio flares is quite peaked while that of the am- 
plitude not. In the framework of the plasmon model 
(! van der Laanl [l966). the width of the flare is related to 
the initial physical size of the plasmon. Thus, the rela- 
tive sharpness of this distribution means that the initial 
radius of the plasmon is roughly always the same. It is 
tempting to interpret that this radius may be constrained 
by the inner radius of the disk. On the other hand, the 
maximum amplitude of the flare is related to the ini- 
tial temperature and energy distribution of the ejected 
plasma. Thus, the relative broadness of the amplitude 
distribution could mean that the amount of energy stored 
inside the plasma prior to the ejection is more variable. 

If we assume that the amount of energy stored inside 
the plasma depends on the duration of the preceding X- 
ray dip, then this relative broadness can be easily ex- 
plained. Indeed, during shorter dips the input of energy 
would be shorter, thus the amplitude of the following 
flare would be smaller. Since the duration of the X-ray 
dips is highly variable, the amplitude of the flares will 
also be variable. On the other hand, the medium which 
stores this energy will always be located in the same re- 
gion, close to the inner radius of the disk; if its physical 
size is constrained by the disk, then the width of the 
flares will be roughly constant. 

One possible interesting model which could provide a 
more precise interpret ation of this beha vior is the "Mag- 
netic Flood" model (T agger et aD2004| ). This model was 
proposed to account for the behavior of GRS 1915-1-105 
during the B class. It relies on the Accretion-Ejection 
Instability (jTagger fc Pellatill999D . which would develop 
during the X-ray dip. This instability relies on the pres- 
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ence of a moderate magnetic field in the inner parts of 
the disk, and can produce Low- Frequency Quasi- Periodic 
OsciUations (LFQPOs), as well as feed a corona with 
matter from the disk. Then, if the magnetic configura- 
tion is favorable, a sudden reconnection event can occur 
between magnetic fields of opposite polarities in the ac- 
cretion disk. This reconnection event would produce the 
spike seen in X-rays, and power the ejection of matter. 
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